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ABSTRACT 

Classical Cepheids are key probes of both stellar astrophysics and cosmology as standard candles and pulsating variable stars. It is 
important to understand Cepheids in unprecedented detail in preparation for upcoming GAIA, JWST and extremely-large telescope 
observations. Cepheid eclipsing binary stars are ideal tools for achieving this goal, however there are currently only three known 
systems. One of those systems, OGLE-LMC-CEP1812, raises new questions about the evolution of classical Cepheids because of 
an apparent age discrepancy between the Cepheid and its red giant companion. We show that the Cepheid component is actually the 
product of a stellar merger of two main sequence stars that has since evolved across the Hertzsprung gap of the HR diagram. This 
post-merger product appears younger than the companion, hence the apparent age discrepancy is resolved. We discuss this idea and 
consequences for understanding Cepheid evolution. 
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1. Introduction 

Classical Cepheids have been crucial for the understanding of 
stellar astrophysics and cosmology since the di scovery of th e 
Cepheid Leavitt law more than a century ago (iLeavitl Il908h. 
They have been used to measure the Hubble constant (iHubblg 
1 192^ iRiess et ^1201 it iFreedman et al.ll?012ll to a precision of 
2% as well as to constr ain the detailed ph ysics of stellar struc¬ 
ture and evolution (e.g.. lBono et al.ll^ObOh . yet there still exist a 
number of mysteries regarding these stars. 

One such challenge is the detailed calibration of the Cepheid 
Leavitt law, i.e., the period-luminosity (PL) relation. In the forth¬ 
coming era of the James Webb Space Telescope, it is expected 
that we will b e able to measure the Hubb le constant to less than 
1 % precision (IFreedman & Madorel2010ll . but this requires mea¬ 
surements of Cepheid fundamental parameters to unprecedented 
accuracy along with independent measurements of Cepheid dis¬ 
tances. The Gaia satellite is currently operating and will measure 
distan ces to thousands of Galactic Cepheids (IWindmark et al.l 
1201 ill , that will complement Large Magellanic Cloud (LMC) 
Cepheids for which distances are known. 

One endu ring uncert ainty is the decades-old ‘Cepheid mass 
discrepancy’ (ICoxlll98^ . where Cepheid masses measured us¬ 
ing stel lar evolution a nd pulsation calcu lations differ by about 10 
- 20% (lKelleill2008l) . lBono et al.l (l2006l) suggested four possible 
resolutions to this discrepancy: changes in the radiative opac¬ 
ities, rotation, convective core ov ershooting in main s equence 
progenitors and stellar mass loss. iNeilson et al.l (1201 ih showed 
that pulsation-d riven mass loss can par tly explain the mass dis¬ 
crepancy, while lAnderson et al.l (l2014ll proposed rotation as an 
alternative solution. Convective core ov ershooting has also bee n 
shown to resolve the mass discrepancy (ICassisi & Salarisll201 111 . 
The solution to the mass discrepancy might simply be a combi¬ 


nation of all three and understanding which physical processes 
are important will constrain both evolution and pulsation mod¬ 
els. However, an ideal method to constrain the Cepheid mass 
discrepancy is to independently measure Cepheid masses. 

Cepheid dis tances and masses can be measured in eclipsing 
binary systems. [Evans et ^ (l200-5ll suggested that about 35% of 
all Galactic Cepheids are in spectroscopic binary systems, but 
none a re known to be in eclipsing binary systems. lEvans et^ 
(1201 3ll presented a catalogue of binary companions detected 
using ultraviolet spectral observations, where the companions 
are all abou t 2 Mp , inferring a larger binary fraction of 60%. 
If vans et^ (l20L5ll report radial velocity measurements which 
further refine the measured binary fraction of Cepheids, noting 
that what has been measur ed is a minimum possible binary frac¬ 
tion. INeilson et al.l (l20L5ll show that a binary fraction of about 
60% is consistent with th e observed binary fraction of ma in se¬ 
quence B-type stars from iKouwenhoyen & de GriisI (l2008l) . 

Th ree eclipsing binary Ce pheids haye been discoyered in the 
LMC dSoszynski et al.ll2008l) . These eclipsing binary systems 
proyide unique mass and distance estimates that can be com¬ 
pared to estimates using Cepheid eyolution and pulsation mod¬ 
els. These help to resolye the Cepheid mass discrepancy and con- 
str ain the PL relation in gre ater detail. 

iPietrzynski et al.l (l2010l) presented detailed obseryations of 
one such eclipsing binary, OGLE-LMC-CEP0227 (CEP0227), 
with a 310 day orbit and a mass ratio between the components 
of g' = 1.00 + 0.01. One binary component is a classical Cepheid 
and the companion a red giant star. The authors measured a 
Cepheid mass M = 4.14 + 0.05 Mp, which i s cons istent with stel¬ 
lar pulsation estimates. ICassisi & SalarisI (1201 ll) compared the 
obseryed fundamental parameters with stellar eyolution calcula¬ 
tions and found that eyolutionary models agree with measured 
parameters if one assumes moderate conyectiye core oyershoot- 
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ing, a result confirmed by other studies (Neilson & Langed2012t 
IPrada Moroni et alJl2012h . IPilecki et~^ (12013 ) analyzed the bi¬ 
nary light curve in greater detail to measure the projection factor 
that is cru cial for using the Baade-Wesselin k method to measure 
distances (lBaadelll92CTlStorm et alJl^l lallSiNgeow et alJ2012t 


iNeilson et al.ll2012cll . This specihc eclipsing binary system con¬ 

strains both the Cepheid mass discrepancy and the calibration of 
the distance scale. 


iGieren et al.l (l2014t) presented new observations of another 
eclipsing binary OGLE-LMC-CEP1718 that appears to be com¬ 
posed of two equal mass Cepheids, both pulsating in the hrst 
overtone. That system has an orbital period of 413-days, but ob¬ 
servations are limited and do not independently constrain each 
star’s radius and mass, making conclusions about the their evo¬ 
lution difficult. This system is appears similar to CEP0227. 

While analysis of CEP0227 helps resolve problems, the bi¬ 
nary system OGLE-LMC-CEP1812 (CEP1812) appears to cre¬ 
ate problems. IPietrzvriski et alj (1201 ih measured the masses and 
radii of the two stars in a 551 day orbit, but found that the 
Cepheid appears to be about 100 Myr younger than its red gi¬ 
ant companion. This result raises questions about the evolution 
of this binary system. While the binary CEP0227 is consistent 
with stellar evolution calculations, the binary CEP1812 is not, 
even th ough the Cepheid mass is consistent with pulsation calcu¬ 
lations. iPietrzvhski et al.l (1201 Ih hypothesized that the Cepheid 
captured the red giant companion into a binary orbit at some 
point during its evolution. 

In this article, we hypothesize that the Cepheid star in the 
systemCEP1812 evolved from the merger of two main sequence 
stars. We discuss our stellar evolution code and models in Sect. 2 
and present measurements of the age difference between the two 
stars in the systems using those models in Sect. 3. In Sect. 4, we 
compute stellar evolution models with changing mass consistent 
with our merger scenario and show that our hypothesis resolves 
the age discrepancy. In Sect. 5 we discuss the implications of a 
stellar merger scenario for understanding Cepheid binaries and 
Cepheid structure and evolution. 


2. Stellar evolution models 


In this work, we com puted stellar evolution models using the 
lYoon & Langed (l2005h code. This code has been used to study 
massive st ar evolution, supernova progenitors an d gamma-ray 
bursts le.g. lCantiello et al.l2009l : lBrott et al.l201 ll) . We have also 
used this code to ex pl ore the evoluti on of classical Cepheids 
dNeilson et alJl2012aH lNeilsonl[2014l) . Models were computed 
assuming moderate convective core overshooting, consistent 
with that m easured for the EMC binary Cepheid OGLE - LMC- 
CEPn227 (ICassisi & Salari^ 120111: iNeilson & Eai^ 120121: 
iPrada Moroni et al.ll2012l) . lit this case, we write the convective 
core overshooting effi ciency as Q.2Hp, that is a fraction of the 
pressure scale height dBrott et alJl201 ll) . Convective core over¬ 
shooting during main sequence evolution acts to create a more 
massive helium core, leading to a more luminous giant star. 
It also acts to prolong main sequence evolution and changes 
the measu red ages of post-main sequence stars. The models in¬ 
cluded the lde .lag er et alJ dl988h ma ss-loss prescription for cool 
stars and the iKudritzki et ak J 19891) recipe for hot stars. We did 
not ass ume enhanced mass l oss during the Cepheid stage of evo¬ 
lution dNeilson et ah 2012al3)- The models a ssume a composi¬ 
tion based on the iGrevesse & Sauvall dl998l) solar abundances 
scaled to the standard EMC metallic ity Z = 0.008 alo ng with an 
initial helium abundance Y - 0.256 dBrott et alJl201 ll) . 
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Fig. 1. Hertzsprung-Russell diagram showing stellar evolution 
models with initial masses = 3.8 M© (grey) and Mi = 
2.7 Mq (black) along with regions consistent with the radius of 
the Cepheid, RcepheW = 17.4 + 0.9 Rq (blue) and that of the red 
giant companion, Rrg = 12.1 + 2.3 Rq (red). 


Stellar evolution models are constrained by the measured 
mass and radius of each component of OGLE-LMC-CEP1812. 
The Cepheid has mass M - 3.74 + 0.06 Mq and radius R - 
17.4 + 0.9 Rq, while the red giant companion has mass Mrq = 
2.64 + 0.4 Mq an d radius Rrg = 12.1 + 2.3 Rq. Unfortunately, 
IPietrzvhski et al.l d201 ll) did not measure the effective tempera¬ 
ture of either star, so we have fewer constrai nts for CEP1812 
than for CEP0227 fe.g. lCassisi & Salarisll201 ll) . 

Given the masses and radii of the two stars in the binary sys¬ 
tem, we computed stellar evolution models with initial masses. 
Ml =3.8 and M 2 = 2.7 Mq. The stellar evolution tracks are 
plotted in Pig. [1] along with the regions of the Hertzsprung- 
Russell diagram consistent with the measured radii. Stellar evo¬ 
lution models appear to ht the measured stellar masses and radii. 
However, the models suggest that the Cepheid has an age of 
about 175 M yr while the red giant star’s age is between 420 
and 450 Myr. fpietrzvhski et al.l (I2OI Ih noted that the stellar ages 
are approximately 190 Myr and 369 Myr for the Cepheid and 
red giant, respec t ively, based on stellar evolution tracks from 
iPietrinferni et al.l (l2004l) . These model age differences are due 
to different amounts of convective core overshooting assumed 
in the models. Including overshooting in stellar evolution mod¬ 
els lengthens the main sequence lifetime and the greater the 
amount of overshooting the longer the main sequence lifetime. 
Regardless of whic h models are cons i dered , there is a signih- 
cant age difference. IPietrzvhski et al.l (1201 ll) suggested the bi¬ 
nary system may have formed by stellar capture. We suggest 
an alternative hypothesis: the apparently younger Cepheid is a 
stellar merger product between two main sequence stars in what 
was originally a triple system that has since evolved across the 
Hertzsprung gap. 

The evolution of the Cepheid in the binary system is also 
notable because the star appears to be crossing the instability 
strip for the hrst time. This phase of evolution is short, about 10^ 
years for stars with mass of about 3-4 Mq. As such, the num¬ 
ber of hrst-crossing Cepheids relative to the number of Cepheids 
evolving on the blue l oop is small, typically about a few percent 
dNeilson et al.ll2012i^ . This fact makes CEP1812 more special 
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as a target for understanding stellar evolution and the transition 
from the main sequence to the red giant branch. 

3. Merger models 


Stellar mergers appear across the HR diagra m, from the 


formation of blue strassler and sdB/O stars dMateo et al. 

199C 

;ISchneider et al.ll2014t). manv massive stars dde Mink et al. 

2013 

), cool R and J stars dizzard et al. 2007; 

Zhane & Jefferv 

2013 

), R Coronae Borealis stars (Clavtonll2013 

) and anomalous 


Cepheids (iBono et al.lll997ll . The coalescence of binary com¬ 
panions can occur via orbital disruptions by third bodies (iKozail 
I l962t iPerets & Kratteij|2012l) o r through tidal interactions and 
Roche Lobe overflow (e.g. lHul]ll98lh . Stellar mergers occurring 
during main sequence evolution quickly settle and evolve in the 
same manner as a main sequence st ar formed with mass s imilar 
to the sum of the merger progenitors dGlebbeek et al.l2()T^ . This 
merger rejuvenates the star, making it appear younger, which 
is consistent with the observed age discrepancy between the 
Cepheid and its red giant companion in CEP1812. 

We propose that the Cepheid component of the eclipsing 
binary system CEP1812 is the result of a merger between two 
smaller-mass main sequence stars that have evolved to become a 
Cepheid. This progenitor system would have been a hierarchical 
triple system in which the red giant was originally the most mas¬ 
sive star. Because of its eccentric orbit (currently about 0.13), the 
outer star (now the red giant) would have induced Kozai oscilla¬ 
tions decreasing the orbital separation between the two smaller- 
mass stars until tides and Roche lobe overflow dominate the evo¬ 
lution. The two stars then coalesced rapidly to form a 3.8 Mq star 
that appears significantly younger than its companion. 

This hypothesis is tested by computing a stellar evolution 
model with mass M\ and adding mass (3.8 M© - M\) after 
about 310 - 330 million years of evolution. At this age, both 
stars are evolving on the main sequence. The mass is accreted 
over a short timescale (relative to the main sequence lifetime 
of about 10 Myr). The main effect of accreting mass is to in¬ 
crease the convective core mass of the star by mixing additional 
hydrogen, hence acts to rejuvenate the star. To first order, this 
mixing is independent of the merging time scale and given the 
masses we consider, addition al mixing effects are not as im¬ 
portant (iGlebbeek et alj|^13h . The merged star then evolves to 
the end of red giant evolution beyond the first-crossing of the 
Cepheid instability strip where Cepheid is observed to be. The 
initial model is assumed to have mass greater than the accreted 
companion, i.e. M\ > 3.8 M© - M\, and the donor and accretor 
have the same surface chemical composition. This is a simplis¬ 
tic calculati on that ignores mass lo st from the system, which is 
0.1 M© dGlebbeek & Polsll2008l) . and ignores issues related 
to mixing caused by accretion. The test is sufficient, however, 
to understand whether a stellar merger resolves the age discrep¬ 
ancy. 

We plot the stellar radius as a function of age in Fig. |2l 
This plot demonstrates the age discrepancy between the 3.8 M© 
Cepheid and its 2.7 M© red giant companion. We also plot stellar 
merger models with progenitor initial masses 2.1, 2.2, 2.3, 2.4, 
and 2.5 M© that each accrete mass until they reach 3.8 M© and 
then evolve normally. Those stellar evolution models are rejuve¬ 
nated and do not cross the Hertzsprung gap until they are much 
older than the main sequence life time of stars with that initial 
mass. The merger models all have a radius consistent with that 
of the Cepheid at approximately the same age as the red giant 
star, implying that CEP1812 is the result of a stellar merger that 
occurred early in the main sequence evolution of a hierarchical 
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Fig. 3. Hertzsprung-Russell diagram showing stellar evolution 
models evolving across the Hertzsprung gap and the first- 
crossing of the Cepheid instability strip. Solid lines follow evo¬ 
lution models with initial mass Mi - 3.8 M© (grey), along with 
stellar merger models with progenitor initial masses 2.1 (violet), 
2.2 (orange), 2.3 (blue), 2.4 (pink) and 2.5 M© (brown). The blue 
region represents the measured Cepheid radius while the dashed 
lines denote the blue edge of the first-overtone instability strip 
‘lOTBE’ (dark blue), the boundary between first-overtone and 
fundamental-mode instability strips ‘FMBE’ (light b lue) the the 
red ed ge of the Cepheid instability strip ‘RE’ (red) (iBono et al.l 

I2000h . 


triple system. However, this solution is not unique, other stellar 
mass combinations are possible depending on when they merge. 

4. Discussion 

There exists a specific class of Cepheids that appear to 
have evol ved from main se quence stellar mergers: anomalous 
Cepheids dBono et al.lll997h . They are typically found in older 
stellar populations in close proximity to RR Lyrae and horizon¬ 
tal bran ch stars, primarily in dwarf spheroidal sa tellite galax¬ 
ies (e.g. iMateo etaPl 199,4 iKinemuchi et al.ll2008h . These stars 
pulsate with periods ranging from about 0.3 to about 2.5 days 
and EMC anomalou s Cepheids have a predicted average mass 
Mac = 1.2 + 0.2 M© dFiorentino & Monellil2012h . Their masses 
tend to be smaller than that of EMC Cepheids and there is a 
small overlap for the longest-period anomalous Cepheid and 
shortest-period c lassical Cepheids. Because of these properties, 
I Sills et aP d2009t) suggested they are the merger product between 
two low-mass main sequence stars that have evolved across the 
Cepheid instab ility strip, akin to the models computed in this 
work, however [Bono et aTI d 19971) showed that a many ACs are 
low-mass core helium-burning stars evolving from the horizon¬ 
tal branch. The anomalous Cepheids do not appear to be related 
to classical Cepheids, but could CEP1812 be a missing link be¬ 
tween these two classes of Cepheids? 

We explore this by plotting our stellar evolution models in 
a Hertzsprung-Russell diagram in Fig. |3 Our merger models 
evolve across the HR diagram with luminosities consistent with 
the 3.8 M© single star model. This implies that CEP1812 is not 
likely to be a traditional classica l Cepheid, but is an anomalous 
Cepheid dFiorentino et alJl2006l) . This possibility is contentious 
because ACs typically have masses about one solar mass and 
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Fig. 2. (Left) Stellar radii as a function of age for stellar evolution models with initial mass Mi = 3.8 (grey) and M 2 = 2.7 Mq 
( black), along with stellar merger models with progenitor initial masses 2.1 (violet), 2.2 (orange), 2.3 (blue), 2.4 (pink) and 2.5 Mq 
( brown). The blue and red colored regions represent the measured Cepheid and red giant companion radii, respectively. (Right) A 
closer view of the stellar radius as a function age consistent with the current age of the red giant 


have lower metallicities. Because of its mass, CEP1812 likely 
has a metallicity similar to the assumed LMC metallicity. The 
star is also not core helium-burning. Both pro perties are incon¬ 
sistent with various pulsation models of ACs (iBono et al.lll997t 
ICaputo et al.l |2004 iMarconi et al.l l2004h . However, a massive 
AC evolving across the Hertzsprung gap would be rare because 
this stage of evolution is short, less than 1 Myr. 

ISoszvnski et akl (|2008|) classified CEP1812 as a classical 
Cephei d based on Eourier d ecomposition of the observed light 
curves dSimon & Leelll98l1) . Measured Eourier parameters al¬ 
low for the classification of variable stars because for differ¬ 
ent variable stars and different pulsation modes the predicted 
Eourier parameters tend to cluster. However, Eourier parameters 
for anomalous Cepheids (ACs) vary significantly and the anoma¬ 
lous Cepheid Eourier parameters from the OGLE-III survey of 
the LMC overlap wit h those of short-period class ical Cepheids 
including CEP1812 (lEiorentino & Monellill2012h . The Eourier 
components for CEP1812 are consistent with both being a short- 
period classical Cepheid or a first-overtone AC. This is some¬ 
what surprising as the Eourier Component (p 2 \ has only a small 
range of values for LMC ACs as measured bv ISoszvhski et aH 
(l2008h . T his result adds credence to the possibility that CEP1812 
is an AC. iPietrzvhski et al.1 (1201 Ih noted that the pulsation pe¬ 
riod and brightness of CEP 1812 is consistent with the measured 
period-luminosity relation (ISoszvnski et al.ll2008h . However, it 
is also c onsistent with the first-overtone AC period-luminosity 
relation dRipepi et al.l2014l) . Both its pulsation period and ampli¬ 
tude are consistent with that of other LMC ACs in the OGLE-III 
survey (ISoszvhski et al.|[2008h again suggesting that CEP1812 
may be an AC and not a classical Cepheid. Based on our mod¬ 
els and the fact that CEP1812 is consistent with the pulsation 
properties of other LMC anomalous Cepheids, we suggest that 
CEP1812 is not a classical Cepheid but is instead an anoma¬ 
lous Cepheid. Thus, CEP1812 is the most massive AC discov- 
ered and is a factor of two more massive than those in the 
lEiorentino & Monellil (1201 2h sample, where the mass is mea¬ 
sured from period-luminosity-color relations. 

It is not surprising that massive anomalous Cepheids are rare. 
Stellar mergers are, themselves, rare, but the most noticeable dif¬ 
ference between a single star evolving and a stellar merger is that 


the latter may appear to be younger. Examples are blue straggler 
stars which are easy to detect in a Globular cluster because the 
majority of stars there are very old. ACs are detected by the same 
method; they tend to be found in old populations along with hor¬ 
izontal branch and RR Lyrae stars. It is this contrast that allows 
anomalous Cepheids to be detected and that also explains why 
they all tend to be about one solar mass. Only low-mass stars 
have life times long enough to both undergo a stellar merger and 
appear in old populations, massive stars that merge would have 
disappeared long before the general population had significantly 
evolved. A merger between a 2.4 and a 1.4 Mq star would not 
appear significantly rejuvenated relative to nearby field stars, but 
does stand out in a binary system. 

Although we suggest that CEP1812 could be classified as an 
AC, it is not clear that the star will have any pr operties signif¬ 
icantly different from other classical Cepheids. iLanged (l2012h 
suggested that strong magnetic fields, typically about 1 kG, that 
are found in ab out 10% of intermediate-ma ss and massive main 
sequence stars dPonati & Landstreeil2009t) could be generated 
in a main sequence merger. If this hypothesis is correct then 
CEP1812 could presently have a weak magnetic field that might 
affect various pulsation properties. Another test of our merger 
scenario is the abundance of CEP1812. If the star is the product 
of a merger of two main sequence stars, where the mass donor 
has underwent some hydrogen burning then one might expect 
some chemical anomalies, particularly in the surface helium, ni¬ 
trogen and carbon abundances. However, because this is a low- 
mass merger, any changes to the carbon and nitrogen abundan ces 
will be small, about a fraction of a dex (iGlebbeek et al.ll^l3h . 

Because the star is a merger product, one might expect some 
abundance anomalies. However, our model assumes that the 
merger is a product of two stars with the same initial compo¬ 
sition, hence any anomalies would be small. While CEP1812 is 
a merger of two main sequence stars, some anomalous Cepheids 
might be created by the merger between a main sequence star 
and an evolved helium-burning star which would generate a 
new helium-burning star with an envelope that is massive rel¬ 
ative to the total stellar mass, or conversely the core mass is 
muc h smaller than for a s ingle star at a similar stage of evolu¬ 
tion (Ide Mink et al.ll2()T^ . These post-main sequence mergers 
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could also evolve to become Cepheids but will have very differ¬ 
ent properties because they have a smaller luminosity for a given 
mass and a more massive envelop leading to different pulsation 
properties. 

In summary, CEP1812 may be an anomalous Cepheid that 
appears to be just like any classical Cepheid such that it may be 
the missing link between the two classes of stars. However, the 
strange evolution of CEP1812 implies that it might be unwise to 
use this star as a calibrator for the Cepheid Leavitt Law or for 
resolving the Cepheid mass discrepancy. 
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